Topological insulator ribbon: Surface states and dynamical response 
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We study theoretically the distributions of charge and spin polarization of a topological insulator 
ribbon, with a realistic rectangular cross section. Due to constriction in two lateral directions, the 
surface states discretize into a series of subbands inside of the bulk band gap. The charge and spin 
distribution show interesting characters which are different from an ideal topological surface state. 
The effect of merging of four different surface states into the new one as an entity are analyzed. 
Optical conductivity and dynamical spin susceptibility of the ribbon are studied. Different from a 
single ideal surface, the optical response and dynamical spin susceptibility of a ribbon do not have a 
clear correspondence. The dynamical spin susceptibility could be used to identify the more adequate 
model for Bi2Se3. 

PACS numbers: 73.20.-r, 78.68.+m, 75.40.Gb 
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I. INTRODUCTION 

The topological insulators are characterized by spin 
polarized helical boundary states protected by time re- 
versal symmetry, [il-fsj In three dimensions, the Bi2Se3 
material has attracted great attention because of its well 
defined Dirac cone like surface state dispersion, 0] 
Though perfect two-dimensional surface states account 
for the most essential properties of topological insula- 
tor, real samples also have lateral pairs of surfaces. For 
strong topological insulators, each surface carries an odd 
number of gapless topological surface modes. :3, 4] When 
the bulk material is constricted in two directions, we get 
a ribbon (or, nano-wire) of topological insulator. |8j] The 
Aharonov-Bohm effect of the Bi2Se3 ribbon showed in- 
teresting (j)Q = h/e period oscillation and a conductance 
maximum at zero fluXjSj, which are in contradiction with 
both diffusive and ballistic transport and are explained 
theoretically in terms of impurity scattering effects. [9l.[l0j 
More conventionally, the ribbon (nano-wire) of Bi2Te3, 
which is also a topological insulator and described by a 
model similar to that of Bi2Se3 6], is considered as a 
promising thermoelectrical material for applications. [llj 
Former theoretical studies of topological insulator rib- 
bon usually start from a cylindrical sample, [ol [lOl [T2I, [Tsj 
While an axi-symmetric configuration is more convenient 
for analytical analysis, a real topological ribbon usually 
has a rectangular cross-section. 8, JLlJ We thus study the 
charge and spin distributions of the surface modes of such 
a realistic topological insulator ribbon. 

For a z • (k X <t) (or equivalently, k • cr) type of ef- 
fective model for the surface states, the charge current 
operator is proportional to the spin operator. The opti- 
cal conductivity for the surface states, which is defined 
in terms of the current-current correlation function, is 
correspondingly directly related to the dynamical spin 
susceptibility. JJJ But for a topological insulator ribbon, 
in which two pairs of lateral surfaces are present, the re- 
lationship between the optical conductivity and the dy- 



namical spin susceptibility would in general be more com- 
plicated. However, the results always depend uniquely 
on the characters of the model used, which also deter- 
mine the properties of the surface states. The study of 
these responses may thus give some criteria to discrimi- 
nate among available models for the topological insulator 
Bi2Se3. 

In a previous work, two different models used in former 
works for Bi2Se3 were pointed out.[l^ Here, we propose 
that the more adequate model for the material could be 
identified by measuring the dynamical spin susceptibil- 
ity of a Bi2Se3 ribbon. While optical conductivities are 
identical, the dynamical spin susceptibilities show clear 
qualitative differences between results obtained from the 
two models. The results are explained in terms of selec- 
tion rule analysis based on direct numerical calculations. 
For optical conductivity, qualitatively different selection 
rules are in effect between transitions along the ribbon 
axis (x direction in this work) and transitions along other 
directions. 



II. MODELS AND METHODS 

Since we are interested mainly in the dynamical re- 
sponses related to the topological surface states, it is 
reasonable to start from the continuum model for a topo- 
logical insulator describing low energy states close to the 
r point of the Brillouin zone (BZ).[6', HI, [i3| 

Formerly, we recognized two different kinds of models 
that exist in the literature for bulk 612X3 (X is Se or Te) 
materials. [ilH-lill Close to the F point, the two models 
could both be written compactly in terms of the Dirac 
matrices as [la 



il(k) -eo(k)/4x4 + ^m,(k)r, 



(1) 



Every unit cell contains two spin and two orbital degrees 
of freedom. The model is hence written in terms of 4 
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by 4 matrices. The two orbitals concentrate mainly on 
the top and bottom (looking along the —z direction) Se 
layer of the various BioSea quintuple units, and are la- 
beled as 1 and 2.[l3,[i3| The basis is taken as V'k=[ciktj 
C2ktj cikj,, C2k^]"^j for a certain wave vector in the three 
dimensional BZ. Since the first term proportional to the 
unit matrix is nonessential to topological properties of 
the system, it is ignored in the following analysis. The 
remaining model is particle hole symmetric. Close to 
the r point, the four remaining coefficient functions are 
mi=o,- ,3(k)={m + \t(kl fc^) ^ tjtl, 3tk^,3tky, 2t^k^}, 
in which t > 0, > and m < (|m| is half of the 
bulk band gap).[l^ In terms of Pauli matrices Si («=0, 
• • • , 3) in the spin subspace and (i=0, • • • , 3) in the 
orbital subspace, the first three Dirac matrices are de- 
fined as@,[ll| TQ=so®ai, Ti=si®a^, T2=S2®o^. Two 
choices of the last Dirac matrix Fa define the two differ- 
ent models in literature, which are: (I)so (X'cr2 and 
(II)S3«)CT3B 20 24]. 

Though a full gap exists in the bulk, gapless modes 
reside on the surface. For surfaces directed along different 
directions, the surface states are described by different 
effective models. Take for example the surface states for 
a sample occupying the lower half z space, the effective 
model for it was solved to beflU 

Heffik.) = M{kxSx + kySy), (2) 

for model I. The charge current matrix is thus j(k) — 
— evk — — eVki?e//(k) — Sets. Where the Pauli ma- 
trices Sx and Sy are in terms of the two basis (?7i)^=(5^i 
and {r]2)i3=Si33. Sap is one for a—l3 and zero otherwise. 
While for model II, the effective model for the surface 
states of the same system isjTsj 

Heff{k) = Mz • (k X s) = 3t{kxSy - kyS^). (3) 

The corresponding charge current is j(k)=— 3eis x z. The 
two bases are ^i = "^[1j 0, 0]^ and 772=-^[0, 0, — i, I]-^. 
For the surface states introduced above residing on the 
xy surface, the two bases have definite spin characters 
for both model I and model II. 

Now consider the surface states of a sample occupying 
the half space y < 0. Similar to surface states on the 
2 = surface, the possible zero energy surface states on 
y — are obtained by solving a set of four coupled second 
order differential equations 

H{kx = 0,ky^ -idy, k, = 0)*(y) = E^{y), (4) 

for E ~ together with the open boundary conditions 
^'(y)|y^o=*(2/)la=-oo=0.[2l|, [111 For model I, the basis 
for the surface states are obtained as 77i = -^[I, 0, 0, 1]^ 

and ?72 = -^[0, 1, — 1, 0]-^ by solving the above coupled dif- 
ferential equations. Now, both orbitals contribute to the 
surface states. The corresponding effective model turns 
out to be 

i/e//(k) = -MkxCTx + "itzk^ay. (5) 



While for model II, the surface states are described by 
an effective model as 

HeffCk) = -'StkxCTx + 2t^k2a^, (6) 

with the same basis set as for model I. The two ba- 
sis states now do not have definite spin characters. On 
the other hand, the Dirac cone of the surface states be- 
comes anisotropic, which is a manifestation of the uni- 
axial anisotropy of bulk states between the xy plane and 
the z direction. fi3\ 

The above effective models for topological surface 
states describe infinite surfaces. A real sample is how- 
ever finite in all three directions. For a ribbon geometry 
realized and studied recently, the sample could be con- 
sidered as infinite along the direction of the ribbon axis 
(taken as x) while finite in the other two directions (taken 
as y and z).^ Then, while k^ could still be taken as a 
good quantum number, real space viewpoint should be 
adopted for y and z dependent quantities. For the sake 
of simplicity, we discretize the y and z coordinates into 
square lattices. The model is thus written as 

3 

Uy .Uz ,kx i—0 

+ J2 {<„Jfcx)[-ro-»r3]i.Vn„«.+i(A:.) + i?.c.} 
+ E {<«.(fc.)[-ro-*r2]^#„,+i,„.(fc.) + H.c.},(7) 

where m-^g ... ^{k)={m + 3t + 2t^ + ^tkl,3tkx, 0, 0}. Uy 
and Uz label the unit cells along the y and z directions. 
The lattice constants along the three directions are taken 
as length units. 

The paramagnetic current operators along three direc- 
tions could be obtained by the continuity equations. [26j 
They are written as 

Jx^-e E V'i„„,(fci:)(3ifca:ro +3tri)'(/'„„„,(fc:r), 
Uy ,nz,kx 

(8a) 

^V^^ E i^tynA^'^)^^y'^riy-l,nAkx) + H.c], {8h) 

jf=e [i;l^^,Akx)iDziPn^^n,^i{kx) + H.C.], {8c) 

where £)j,=-|<(ro -ir2) and £'z=-tz(ro -iFa). In this 
work, we would focus on the zero temperature response 
behaviors of the system. The dynamical conductivity is 
obtained from the retarded current-current correlation 
functions by the Kubo's formula. [25!, l2^ Ignoring contri- 
bution from the diamagnetic current (which contributes 
to the zero frequency Drude weight), the expression is 

o-a/3(a;) = — nQ^(w), (9) 
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where the retarded current-current correlation function 
is defined as 

• I-+00 

n„M^) = / dte-*0W(o|[j-rw,j^(o)]|o) 

^ J — oc 

V + Eo-En+iv u; + En - Eo + iv ' 

{\n)} is a complete set of basis states, with |0) denoting 
the ground state. En is energy of the state \n). rj is 
the positive infinitesimal and taken as a small positive 
number in realistic calculations. 

The spin operator of the system is defined as S 
= k T^nyu, J2k^ ^lyn, {kx)&®(To^nyn, {K)- The zcro tem- 
perature dynamical spin susceptibility (for q = 0) is thus 
defined as 

x.,H = -- j ^ dte:-'e{tm[s,{t),s,{Q)m 

_ 1 (0|5,|n)(n|^,-|0) (0|5,-|n)(n|^,|0) 

V ^^Lj + Eo-En+iri u + En-Ea + iif^ ' 



III. RESULTS AND DISCUSSIONS 

Without loss of generality, we consider a ribbon with 
Ny^N^=2Q and for t==t^==0.5, to=-0.7. Fig. 1(a) is the 
dispersion of the ribbon. Within the bulk gap region 
(|£'| < —m), sub-bands form with a finite minimal gap as 
a result of discretization of the topological surface states 
induced by lateral constriction. As shown in Fig. 1(a), 
the sub-bands are labeled symmetrically with respect to 

= as {±1, ±2, . . .}, with positive (negative) sign 
denoting positive (negative) energy sub-bands and the 
larger numbering representing a larger distance to i? = 
0. Every state is twofold degenerate due to combined 
time reversal symmetry and inversion symmetry of the 
material. 

The charge and spin distributions of the topological 
surface states are essential to its physical properties. In 
the ribbon geometry, an interesting feature is related to 
the hybridization and redistribution of different surface 
states. We give the spin and charge distributions for sev- 
eral typical surface modes on the sub-bands within the 
bulk gap. As shown in Fig. 1(b) for the charge distri- 
bution of a typical surface mode for kx=0, the charge is 
distributed centro-symmetrically close to the boundary 
of the ribbon's cross-section. Around the cross-section, 
more charge is distributed close to the four corners which 
is easily understood from the smaller effective radius 
there compared to the flat part of the sample. As shown 
in Figs. 1(c) and 1(d), the charge distributions for the 
twofold degenerate modes of n—1 for fc^^ClTr are asym- 
metric and centered around one corner of the ribbon, 
displaying clearly the effect of hybridization of two kinds 
of surface modes across the corner. 

While the charge distributions are qualitatively identi- 
cal along the y and z directions and are the same for the 
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FIG. 1: (a) Dispersion of a topological insulator ribbon with 
Ny=N^=20. (b) Charge distribution (multiplied by 10^) of 
the surface modes n=3 for kx~0, results for other n (for kx=0) 
within the bulk gap region are similar, (c) and (d) are the 
charge distribution (multiplied by 10^) of the two fold degen- 
erate n=l surface eigenmode for kx—O.l-K. 
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FIG. 2: Spin distribution for a surface mode labeled as the 
first sub-band at kx—O. The three components of the spin 
expectation value on every site are presented in (a), (b) and 
(c) (d, e and f) for model I (model II). The results in (a), (b) 
and (c) (in d, e and f) are multiplied by 10"^ (10*). 
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FIG. 3: Longitudinal dynamic conductivity (optical conduc- 
tivity) along X, y and z directions for three typical chemical 
potentials (a) fj. = 0, (b) fi = 0.3, and (c) fi = 0.7. Results 
for the two models are identical. 2jm| = 1.4 is the bulk band 
gap. 77=0.01. 



two different models, the spin distribution is more com- 
plicated. For a state residing on an infinite xy surface, 
the spin only has in-plane component. While for states 
on an infinite xz surface, the basis of the surface states 
indicate that the spin has both in plane and out of plane 
components. In the present ribbon geometry, the surface 
modes resulting from hybridization of the two types of 
surface modes would be different from an ideal surface. 
Fig. 2 shows distributions of the three spin components 
for a typical surface mode, which is one of the two degen- 
erate modes for n=l and kx—O- Figs. 2(a, b, c) are for 
model I. A peculiar feature is that the expectation value 
of the y spin component is everywhere zero. However, 
spin distributions of the corresponding surface state for 
model II shown in Figs. 2(d, e, f) have all the three spin 
components nonzero along the edge of the cross-section. 

Having the above charge and spin distributions of the 
surface sub-bands in mind, it is natural to expect that 
the charge response of the two models would be similar, 
but the spin responses would possibly show very differ- 
ent behaviors. Fig. 3 shows the optical conductivity (in 



unit of e'^/h) of a ribbon with N^=20 and Ny=20. The 
same results are found for the two models. The optical 
conductivity is defined as real part of the longitudinal 
dynamical conductivity and is related to optical absorp- 
tion. One interesting feature is the low energy peaks 
of optical absorptions related to optical transitions be- 
tween the surface sub-bands inside of the bulk band gap 
for oj < 2\m\ = 1.4. [23| When the chemical potential lies 
at the charge neutral point, only the x component of the 
optical conductivity is significant. The y and z compo- 
nents of the optical conductivity are smaller and nearly 
featureless. When the chemical potential is shifted away 
from the Dirac point by A, the low frequency part be- 
low 2|A| for <Txx is depleted while the higher frequency 
part remains unchanged. In contrast, some low frequency 
(a;<^2|m|) peaks emerges for ayy and azz- This reflects 
a difference in the selection rules of optical transitions 
along the three directions. Since is still a good quan- 
tum number, only vertical optical transitions between 
states with the same kx are allowed along all three direc- 
tions, as could be seen from the current operators defined 
in Eq. (8). 

Differences between axx and optical conductivities 
along the other two directions arise from selection rules 
with respect to n in the optical transitions. Numerically, 
only transitions between pairs of sub-bands labeled by 
n and — n are allowed for optical transitions along the 
X direction. While for optical transitions along y and z, 
the transitions involving sub-band n is largest with re- 
spect to n ± 1. So, while for axx a series of peaks below 
2\m\ are expected corresponding to transitions between n 
and —n sub-bands, (Tyy and azz only have low frequency 
peaks corresponding to transitions between adjacent sub- 
bands. For nonzero A, the low frequency (w < 2|A|) op- 
tical transitions contributing to strictly forbidden 
since the involved sub-band pairs are all occupied. Thus 
the corresponding low frequency part of axx is depleted. 
But for ayy and azz , the number of states that could con- 
tribute to low frequency optical absorptions continuously 
increases, resulting in a monotonic enhancement of the 
low frequency (w < 2|A|) peaks with increasing |A|. This 
nicely explains the evolution of the low frequency optical 
conductivities in Fig. 3. 

We now proceed to calculate the dynamical spin sus- 
ceptibilities. In accordance with the optical conductivity, 
we focus on the imaginary part of the dynamical spin sus- 
ceptibilities. Only the longitudinal components would be 
analyzed here. For an ideal xy surface, the optical con- 
ductivity is identical to the dynamical spin susceptibility. 
In the presence of xz surface, another correspondence be- 
tween optical conductivity and spin susceptibility specific 
to the xz surface states mixes in. In Fig. 4 we show the 
imaginary parts of the longitudinal dynamical spin sus- 
ceptibilities normalized by the frequency. In agreement 
with expectation based on the different spin distributions 
of the surface sub-bands, the dynamical spin responses 
for the two different models show qualitatively different 
behaviors and thus could be used as a means to tell the 
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FIG. 4: Dynamic spin susceptibility, the imaginary part nor- 
malized by the frequency, (a,) fi = for model I, (b) fi — 0.3 
for model I, (c) = for model II, (d) fj, = 0.3 for model II. 
2\m\ — 1.4 is the bulk band gap. ry=0.01. 

right model for Bi2Sc3. 

Though the two types of surfaces are coupled together, 
some qualitative correspondences between the optical 
conductivity and the dynamical spin susceptibility could 
still be identified, [l^ In particular, for an infinite xy sur- 
face, the X component of the optical conductivity is pro- 
portional to the X component of the dynamical spin sus- 
ceptibility for model I while it is proportional to the y 
component of the dynamical spin susceptibility for model 
II. Qualitatively similar correspondence could be identi- 
fied from comparing Fig. 4 with Fig. 3. The y and 
z components of the optical conductivity suffer stronger 
influences from the lateral constriction and the induced 
surface states hybridization, the correspondences to the 
dynamical spin susceptibilities are also not very straight- 
forward. Note that, despite the similarity with the x 
component of the optical conductivity, a subtle difference 
is that the first peak of the dynamical spin susceptibility 
occurs at a frequency which is twice that of the first peak 
in axxii-^)- This reflects a general difference between se- 
lection rules for electronic dipolar transitions (contribut- 
ing to (TctQ, (cli)) and magnetic dipolar transitions (con- 
tributing to Xaa{uj))^ which coincides only very rarely 
for a state such as the surface state on the xy surface. 
Explicitly, the matrix elements for the spin operators are 
nonzero only for transitions connecting second neighbor 
sub-bands, that is between n with n ± 2, profoundly dif- 
ferent from the above-mentioned selection rules for the 
optical conductivity. 

We would like to point out an intricate point on iden- 
tifying the right microscopic model for Bi2Se3. Exper- 



imentally, the spin polarization of the xy surface states 
was found to be perpendicular to the 2D wave vector, psj 
Though model II gives naturally the right spin polariza- 
tion, it could not directly be chosen as the right model. 
Changing the spin-orbital coupling in the xy plane from 
k ■ (T to z • (k X cr) in the bulk model, the same spin 
polarization as observed in experiment is realized also 
for model I. Since this substitution amounts to a 

simple redefinition of the spin axes in the xy plane, Xxx 
and Xyy simply exchange with each other and thus simi- 
lar distinction between the two models by the dynamical 
spin susceptibilities is still feasible. 

Finally, we mention the experimental measurements of 
the dynamical spin susceptibilities. The electron mag- 
netic resonance (EMR) could approach ~3 THz which 
amounts to approximately 12 meV. [291] For a ribbon with 
much larger cross-section than considered above but sim- 
ilar to experiment, [s^ the sub-bands inside of the bulk 
band gap would be much denser than for Ny = = 20. 
The corresponding low frequency peaks in the dynam- 
ical spin susceptibility would also within the reach of 
EMR. [29] Another technique is the inelastic neutron scat- 
tering (INS). Though it is usually hard to see surface 
effects in terms of INS, when the chemical potential is 
tuned close to the Dirac point inside of the bulk band 
gap, the surface would possibly give the dominant sig- 
nal. At last, the spin flip surface Raman scattering could 
be another possible way to measure the dynamical spin 
susceptibility. [sOl 



IV. SUMMARY 

We have calculated the charge and spin distributions 
and the corresponding dynamical spin and charge re- 
sponses of a topological insulator ribbon. Two models 
giving different surface states are analyzed as a compar- 
ison. Constriction along two lateral directions change 
the gapless Dirac cone like surface states into sub-bands 
inside of the bulk band gap. The corresponding charge 
distributions of the sub-band states resulting from hy- 
bridization of four different surface states are identical for 
the two models, giving rise to identical optical conductiv- 
ities for the two different models. The spin distributions 
of the surface states are however quite different between 
the two models. The dynamical spin susceptibilities are 
thus also quite different and could be used to identify the 
right model for Bi2Se3. 
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